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. Abstract. At present three pulsars are known which 
clearly show a strongly increasing degree of circular po- 

. larization with frequency. As this is accompanied by a 
smooth decrease of linear polarization, we investigate, if 

I this observation can be explained through a propagation 
mode transition from linear to circular polarization. Us- 
ing a previously published model we find that the rate 
of change of the polarization types with frequency is well 

. consistent with the theory. 

I The small number of only three objects does of course 

■ not allow to draw a general conclusion for the whole sam- 
ple of pulsars. However we show that the very unusual 

. behavior of these three objects can be well explained with 
this model and we discuss, why only so few such objects 

. are presently known. 
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1. Introduction 

Whether the observed polarimetric features of pulsar ra- 
dio emission are intrinsic properties of the emission mech- 
anism or if they are a result of propagational effects is one 
of the open questions for the understanding of the polar- 
ization of radio pulsars. Early observations of the 'S'-type 
swing of the polarization position angle through the pulse 
profile of the Vela pulsar led to an geometrical interpre- 
tation by Radhakrishnan & Cooke (1969). In their model 
the swing reflects the geometry of the magnetic field lines 
projected on the plane of the sky. In some cases the ex- 
pected swing is matched so well that a non-geometrical 
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origin for the position angle seems rather improbable. Ru- 
derman & Sutherland (1975) gave an explanation for this 
phenomenon as an intrinsic property of emission by co- 
herent curvature radiation. The particles experience an 
acceleration perpendicular to their motion in the plane of 
curvature of the field line. The resulting radiation has its 
electric field vector within this plane, thus producing the 
geometrical signature. 

Alternatively, it was proposed that the observed po- 
larization originates from propagational effects or is at 
least influenced by them. Various authors calculated the 
properties of the propagational modes in a pulsar mag- 
netosphere (e.g. Allan & Melrose 1982, Barnard & Arons 
1986, Lyutikov 1998 and references therein). It is widely 
agreed among those authors, that two independent, gen- 
erally elliptical polarization modes propagate, which are 
oriented parallel and perpendicular to the plane of curva- 
ture of the magnetic field line. The observed polarization 
then corresponds to the shape of these modes at the dis- 
tance from the neutron star, where the radiation decouples 
from the plasma. This distance is called the polarization 
limiting radius (hereafter PLR) and does not necessarily 
coincide with the place of emission. The magnitude of the 
PLR has been a subject of theoretical debates. Barnard 
(1986) and Beskin et al. (1993) place the PLR at the cy- 
clotron resonance. Melrose (1979) defines a coupling ra- 
tio Q — l/(AfcL) where Afc denotes the difference in 
wavenumber of the two modes and L is a characteristic 
length scale. The PLR is at the position where Q ~ 1. 

In a previous paper (von Hoensbroech et al. 1998b) 
we have shown that many of the complex variety of radio 
pulsar polarization states can be understood qualitatively 
if some propagational effects in the pulsar magnetosphere 
are taken into account. This was done by using a simple 
approximation for the properties of natural polarization 
modes as they propagate through the magnetosphere. The 
model is based on the following assumptions: 1. The tem- 
perature of the background plasma is assumed to be zero 
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(distribution function 7(7) = (5(7 — 7bg), hereafter 7bg is 
the 7-factor of the background plasma). 2. Furthermore no 
strong pair production in a sense, that the e~e"'"-plasma 
is not approximately neutral, is assumed. Considering all 
relevant angles under which the propagating wave and the 
magnetic field lines intersect, the shape of the natural po- 
larization modes can be calculated at any given point be- 
tween the emission height and the light cylinder radius 
(hereafter Rlg, Rlc — cP /2tt). 3. Assuming a PLR and 
a certain Lorentz-factor for the background plasma, the 
following qualitative statements can be made: a. The po- 
larization modes propagate independently with their po- 
sition angle parallel and orthogonal to the local plane of 
field line curvature, h. the radiation depolarizes towards 
high frequencies, c. the degree of polarization correlates 
with the pulsars loss of rotational energy E and d. the po- 
larization smoothly changes from linear to circular with 
increasing frequency. 

This last point was met by recent observations at the 
relatively high radio frequencies of z/obs = 4.9 GHz (von 
Hoensbroech et al. 1998a) as it was already suspected by 
early measurements of Morris et al. (1981). For a few pul- 
sars the degree of circular polarization increases strongly 
with frequency. For three of those objects, the circular po- 
larization reaches values of more than 50 %, even surpass- 
ing the linear polarization. This rather unusual behavior 
appears smooth and is thus likely to be systematic, as an 
inspection of lower frequency data showed. We therefore 
use the data of those three clear examples to quantita- 
tively compare the rate of change from linear to circular 
to the predicted one in statement d. 

2. Analysis and Results 

Presently we have data for three pulsars, where a clear 
change from linear to circular polarization towards high 
frequencies is observed (see Fig. 5 in von Hoensbroech et 
al. 1998b). It is obvious from Fig. 1 that the three ob- 
jects have very different rotational periods P and period 
time derivatives P. Hence they do not form an isolated 
group with respect to their rotational parameters. Since 
the change from linear to circular polarization is super- 
posed by general depolarization effects, which affect both 
types of polarization likewise, the degrees of polarization 
cannot be compared directly with the theory. Therefore 
we choose the ratio 



10' 



10' 



I 10" 



0) 
Q 
■a 

% 10 

0) 
Q. 



10' 



10^ 



B1 737-30 



B1913+1 



B01 44+59 




10^ 



10' 



10' 
Period [s] 



10" 



10' 



Fig. 1. Position of the three pulsars in the P — P-diagram for 
which the change from linear to circular polarization has been 
observed. They do not form a coherent group within the pulsar 
sample. One of them is very young, one is rather average and 
the third one is a relatively old object. 



Table 1. Polarization data and references for all three pulsars. 



Pulsar 


Freq. [MHz] 


nL[%] 


nc[%] 


R 


Rcf.tl 


B0144+59 


610 


29 ±2 


11 ±3 


2.64 ±0.74 


[1] 




1408 


23 ± 1 


24 ±2 


0.96 ± 0.09 


[1] 




1642 


25 ±6 


26 ±6 


0.96 ±0.32 


[1] 




2695 


16 ±3 


50 ±4 


0.32 ±0.07 


[4] 




4850 


12 ± 1 


50 ±2 


0.24 ±0.02 


[2] 


B1737-30 


610 


22 ± 1 


20 ± 1 


1.10 ±0.07 


[1] 




660 


43 ±3 


17 ±4 


2.53 ±0.62 


[3] 




925 


73 ±2 


21 ± 2 


3.48 ± 0.34 


[1] 




1408 


80 ± 1 


38 ± 1 


2.11 ±0.06 


[1] 




1642 


71 ± 1 


45 ± 1 


1.58 ±0.04 


[1] 




4750 


53 ±7 


59 ±5 


0.90 ±0.14 


[2] 


B1913-hl0 


408 


18 ±2 


13 ±3 


1.38 ±0.35 


[1] 




610 


38 ±2 


22 ± 2 


1.73 ±0.18 


[1] 




925 


77 ±4 


35 ±3 


2.20 ±0.22 


[1] 




1408 


50 ± 1 


42 ± 1 


1.19 ±0.04 


[1] 




1642 


49 ±2 


46 ±3 


1.07 ±0.08 


[1] 




4850 


21 ±4 


55 ±5 


0.38 ±0.08 


[2] 



R := 



Hl 
He 



(1) 



between the degrees of linear (Hl) and the circular (He) 
polarization as an depolarization independent parameter. 
R and its statistical error Ai? can easily be calculated 
from the data. 



^ Reference-code: [1] Gould & Lyne (1997), [2] von Hoens- 
broech et al. (1998a), [3] Qiao et al. (1995), [4] unpublished 
Effelsberg data. 



The polarization data we used for this analy- 
sis were accessed through the online EPN-database 
( fittp://'ww'w.mpifr-bonn.mpg. de/pulsar/ data/ ) . The cali- 
brated data were extracted in the EPN-format (Lorimer 
et al. 1998). The degrees of polarization were calculated 
using the same routine for all profiles. All relevant param- 
eters and references are listed in Table 1. 

The theoretical functional dependence of R can be 
derived from equation (25) in (von Hoensbroech et al. 
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1998b). Applying a Taylor-expansion in inverse frequency 
to the first order yields 



R{u) 



,PLR 



sin^ 6'pLR 



7bg 
+0 {i 



47r(cos 

-3^ 



^plr-/3)(1-/3cos0plr)2 



(3) 



Here ile.PLR is the local electron gyro frequency at 
the PLR, 0pLR the angle between the propagating wave 
and the direction of the local magnetic field at the PLR 
and (3 = v/c = y^l — l/T^g the streaming velocity of the 
plasma. Higher order terms can be neglected as long as 
the wave frequency is different from ~ ile/STT. 

Obviously, a -frequency dependence of R is re- 
quired by the theory. A comparison of the observed fre- 
quency dependence of R with the predicted one is there- 
fore a strong test for the theory. Figs. 2-4 show a compar- 
ison of the data points and the theoretical curve (dashed 
line) for R. 

There are a couple of parameters which indirectly en- 
ter Eq. (2) as scaling factors. fie.PLR is proportional to the 
local magnetic field strength B{PLR). This value again 
depends on basic pulsar parameters such as the period, 
its time derivative and, if known, on the inclination angle 
between the rotation- and the magnetic dipole axis. Fur- 
thermore the angle ^plr between the propagating wave 
and the local magnetic field depends on the chosen field 
line and the assumed emission height. For the background 
plasma Lorentz factor we made the assumption 7bg = 200. 
Finally we chose the PLR at 20 % of i?LC- As the values 
of the PLR and the emission height (2 % of i?LC are given 
as fractions of i?LC, their absolute values depend on the 
period. The combination of 7bg and PLR was chosen with- 
out restriction of generality as various other combinations 
yield the same result (see Fig. 5). 

However, apart from the known intrinsic pulsar param- 
eters P and P, the same set of parameters was used for all 
three pulsars. Please note that these parameters only enter 
Eq. (2) as scaling factors, yielding a parallel displacement 
of the function. Hence they determine the frequency range 
where the transition from linear to circular polarization 
takes place, but not the functional dependence ~ v~^. 

2.1. PSR B0144+59 

This pulsar is the first one in which we found the effect 
of increasing circular polarization to high frequencies. Its 
rotational parameters are P = 196 ms and P = 0.3 x 
10""^^, yielding to a weak surface magnetic field of 'only' 
Bo = 2.3 X 10^ T, an average value for ^ = 1.3 x lO^^ W 
and the characteristic age r = 1.2 x lO'' yrs. Reasonable 
data in full polarization was available between 610 MHz 
and 4.85 GHz. 

Fig. 2 shows the measured values and the theoretical 
curve for the change of R with frequency. 
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Fig. 2. PSR B0144+59, R = Ul/Uc versus frequency. The 
theoretical change of R (dashed line) is compared to measured 
data. The horizontal dotted lino corresponds to i? = 1, i.e. 
linear and circular polarization are of equal strength (here 
R = 1 at v = 1.4 GHz). Parameters for theoretical line: 
7i,g = 200, PLR = 20% of Elc. 
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Fig. 3. PSR B1737-30, sec caption of Fig. 2 and text for do- 
tails (here R = 1 at v = 3.5 GHz). See text about the two "out- 
riders'. Parameters for theoretical line: 7bg = 200, PLR = 20% 

of i?LC. 



2.2. PSR B 1737- 30 

PSR B 1737-30 has a spin period of P = 607 ms and a 
period derivative of P = 466 x 10~^^. The resulting spin 
down energy loss E = 8.3 x lO^^H^ places it amongst the 
top 10% of the pulsar sample. The very low characteristic 
age r = 2 X 10"* yrs and the very high surface magnetic field 
Bo = 1-7 X 10^ T make this one an extreme object. Note 
that in terms of the surface magnetic field, this object is at 
the opposite end of the "normal" pulsar sample compared 
to PSR B0144+59. 

The change of R with frequency is shown in Fig. 3. 
Although the two low frequency points do not fit the the- 
oretical function, the frequency change of R smoothly fol- 
lows this function at higher frequencies. Note that the 
two outrider profiles are significantly affected by interstel- 
lar scattering which could also alter the polarization state 
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Fig. 4. PSR B1913+10, see caption of Fig. 2 and text for 
details (here R — 1 ai u = 1.8 GHz). Parameters for theoretical 
line: 7bg = 200, PLR = 20% 



2.4- Outriders 

The systematic deviation of the low frequency points is 
certainly a draw back of these observations. However they 
can be understood through the following argument: von 
Hoensbroech et al. (1998a) have shown that the polariza- 
tion properties in general are much less systematic at low 
radio frequencies compared to higher ones. This indicates 
that the polarization of pulsars undergo some sort of ran- 
domization at low radio frequencies. This can be caused 
either through intrinsic variations - e.g. non constant PLR 
at low frequencies - or through additional propagation ef- 
fects in the highly magnetized medium close to the pulsar, 
which mainly affect low radio frequencies. 



3. Discussion 
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Fig. 5. Possible combinations of PLR and 7bg which yield 
consistent results for the prefactor in Eq. (2) which fit to the 
observed data of each pulsar. A variation of this factor results 
into a parallel displacement of the theoretical curves in Figs. 
2-4. The solid line below corresponds to the assumed emission 
height at 2 % of the i?LC as the minimum PLR. 



when emission from different pulse phases with different 



polarization states are superposed (see also Sect. 2.4). 



2.3. PSR B 191 3+ 10 

With a spin period P — 405 ms and its temporal derivative 
P = 15 X 10""'^^ this is an average pulsar. This is also 
reflected by its parameters E = 9x 10^'^ W, Bq = 2.5 x 10^ 
T and r = 4 x 10^ yrs. 

Fig. 4 shows the measured values for R and the theoret- 
ical function . As for PSR B1737— 30 the lower frequency 
points are too small, but at higher frequencies the theoret- 
ical curve is matched perfectly. As for the previous pulsar, 
the two outrider profiles are significantly scattered. 



Based on a simple propagational scenario for natural po- 
larization modes in pulsar magnetospheres we were able 
to show that for the available data set of three pulsars, 
theory and observations are in very good agreement. In 
these objects the polarization changes smoothly from lin- 
ear to circular with increasing frequency. The principal 
theoretical i/~^-functional dependence of of this change, 
expressed in Eq. (2) is independent of all pulsar param- 
eters (as long as a mono-energetic background plasma is 
concerned). The individual pulsar parameters enter Eq. 
(2) only as scaling factors, hence yielding a parallel dis- 
placement of the function. Thus the comparison of the 
calculated and observed functional dependence is a strong 
test for the theory. As shown in Sect. 2 the rate of change 
of R with frequency is consistent with Eq. (2). 

We note that for all three pulsars the same combina- 
tions of 7bg and PLR result into quantitative agreements 
between the measured points and the theoretical curve, 
although these objects are very different in their param- 
eters P and P. Obviously not only the dependence 
is in agreement with the observation, but also the scaling 
of the function as it is defined through the prefactor in 
Eq. (2) for these three pulsars. This scaling coincidence is 
remarkable but does not necessarily mean, that all pulsars 
have to scale in the same manner. 

Of course the number of pulsars which clearly show 
this effect is too small to make a statistically rigorous 
statement. This is especially true as it cannot be pre- 
dicted at which frequency this effect should occur and 
hence for which other pulsars we should have observed 
it. We emphasize that it is rather difficult to observe such 
a frequency behavior of polarization since it requires a 
significant increase of the circular polarization to result 
into some decrease of the linear because they are added in 
quadrature. 

Pulsars exhibit steep radio spectra, a fact which re- 
duces the total number of available profiles in full polariza- 
tion to about 100 at 4.9 GHz with many of them at a fairly 
low signal-to-noise ratio. As pulsars are also known to de- 
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polarize towards high frequencies, the signal-to-noise-ratio 

of the polarized emission is often very low, hardly allowing 
a credible determination of the ratio between linear and 
circular polarization. Nevertheless, it has been shown at 
low frequencies that the linear polarization reduces with 
frequency while such a trend is not evident for circular po- 
larization (e.g. Han et al. 1998). This result qualitatively 
fits into the same trend of changing the ratio between the 
two polarizations. Taking all points mentioned above into 
account we are not surprised that only three pulsars have 
been found yet which clearly show this eflfect and allow an 
opportunity to test the theory. 

In this paper we have used the frequency dependence 
of integrated profiles. These integrated profiles of course 
differ strongly from the individual pulses. It has been 
shown e.g. by Stinebring et al. (1984), that these indi- 
vidual pulses are often highly circularly polarized also at 
low frequencies. This is not necessarily in contradiction to 
the ideas described in this paper as the conditions within 
the magnctosphere arc known to be highly non-stationary 
and unstable. The location at which the radiation decou- 
ples from the plasma (PLR) might therefore vary strongly 
from pulse to pulse. However, the integrated profile repre- 
sents an average over these individual pulses, correspond- 
ing to an average PLR. A crucial test whether the de- 
scribed effect is responsible for the circular polarization 
would be the simultaneous observation of single pulses in 
full polarization over large frequency intervals. 

Summarizing, we may say that the observations of the 
three objects support the theoretical approach that the 
polarization characteristics - and especially the high fre- 
quency occurrence of strong circular polarization of at 
least these three pulsars are explainable in terms of prop- 
agating natural wave modes in the magnctosphere. Thus, 
it may be worth to follow that line of thought from the 
theoretical point of view and to try to get some more high 
quality data of the general polarization properties at low 
and high radio frequencies. 
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on the manuscript and support of this work. 
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